Rigorous Proof of Pseudospin Ferromagnetism in Two-Component Bosonic Systems 

with Component-Independent Interactions 
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For a two-component bosonic system, the components can be mapped onto a pseudo-spin degree 
of freedom with spin quantum number S — 1/2. We provide a rigorous proof that for a wide- 
range of real Hamiltonians with component independent mass and interaction, the ground state is 
a ferromagnetic state with pseudospin fully polarized. The spin-wave excitations are studied and 
found to have quadratic dispersion relations at long wave length. 
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There has been an explosion of interest in the physics 
of interacting bosonic systems since the realization of 
Bose-Einstein condensation in trapped alkali gases 0. 
Recently, much attention has focused on bosons with in- 
ternal degrees of freedom. These discrete internal degrees 
of freedom come from the (nearly) degenerate hyperfine 
levels or other sources that give rise to multiple compo- 
nents for the bosons. Indeed, recent studies [H,||J|,||,^,0,|| 
demonstrated that the multi-component bosonic systems 
have much richer physics than their single-component 
counterpart due to the additional degrees of freedom. 

Despite intensive research, there has been relatively 
few rigorous results on these interesting systems. In this 
paper we study an example of such systems, namely 
a two-component bosonic system with component- 
independent mass and interactions Q . Obviously the two 
components of the boson can be mapped to a pseudo-spin 
degree of freedom with spin quantum number S — 1/2 
(i.e., pseudo-spin "up" state representing one compo- 
nent, pseudo-spin "down" state representing the other 
component), and the interaction is pseudo-spin indepen- 
dent, thus the system is invariant under pseudo-spin 
SU{2) rotation. We prove rigorously that for a large class 
of real Hamiltonians, the ground state of the system is 
a fully polarized pseudo-spin ferromagnet. The origin of 
the ferromagnetism is neither spin-dependent interaction 
(as there is none in our model), nor interparticle inter- 
action as Coulomb interaction driven ferromagnetism in 
electronic systems; instead it is due to the fact the ki- 
netic energy term of a bosonic system forbids the ground 
state wave function to change sign, thus forcing a totally 
symmetric spatial wave function and hence a totally sym- 
metric (and fully polarized) pseudospin wave function. 
Thus the ferromagnetism is driven by kinetic energy in 
our system. We will also study the ferromagnetic spin 
wave spectra of the system under various conditions. For 
simplicity we will refer to the pseudospin of the bosons 
as spin from now on. 

Consider the following Hamiltonian describing N iden- 
tical bosons: 
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Here U is the trapping potential and V is the two- 
particle interaction. Our proof can also be generalized 
to cases with multi-body (three and above) interactions 
when present. Since H is spin independent, we have 
[H, Si] = 0, where is the spin operator of ith particle. 
In particular, it possesses the global SU{2) symmetry: 
[H, Stot] = 0, with Stot — J2i Si- Thus the eigenkets of 
H may be chosen to be simultaneous eigenkets of S^^j 
and = (TVt - Ni)/2. 

We now proceed by considering the ground state of 
an enlarged Hilbert space for H, namely that of N dis- 
tinguishable particles, i.e., no permutation symmetry is 
required for the wave function. Since H is spin indepen- 
dent and commutes with S^ for every 1 < i < iV, in this 
enlarged Hilbert space we can always choose the eigen- 
kets of H to be simultaneous eigenkets of Sf , so that the 
wave functions take a factorized form llOtl: 
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Since the eigen energy is independent of the spin wave 
function, wc only need to focus on the spatial wave 
function of the ground state: </)o(ri, r2, • • • , rjv). Since 
H is real, (j>Q can be chosen to be real as well. We 
now prove that (pQ can be chosen to be real and 
non-negative, using a well-known trick. If (pQ con- 
tains both positive and negative parts, we can then 
construct a trial wave function that is non-negative: 
0g(ri,r2, • • ■ ,rAr) = |0o(ri,r2, • • ■ ,i"Af)L and show that 
{(l)o\H\(j)o) = {(j)o\H\(po) . It is easy to see that the expecta- 
tion values of the potential terms are the same; the expec- 
tation values of the kinetic energy term are also the same 
because the integrand is the exactly the same everywhere 
(f>o 7^ 0, while the singularity of V^^o where 0o is not 
strong enough to overcome the speed that cpo and 
makes no contribution. Thus we find from variational 
principle that the ground state energy Eq < {(f>o\H\(j)o) . 
The equal sign holds only when (pQ happens to be the 
true ground state, which is not the case under generic 
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situations. Thus generically the ground state wave func- 
tion 00 is non-negative definite, and moreover, positive 
definite for non-singular U and V . This also implies 
that the ground state must be non-degenerate, because 
if there existed two different ground state wave functions 
that are both non-negative, it is always possible to con- 
struct new ground state wave functions that have both 
positive and negative parts by making linear combination 
of the two. 

Having established that 0o is non-degenerate and non- 
negative definite under generic situations, it follows that 
it must also be totally symmetric: Pij4>o — (t>Oi where Pij 
is the permutation operator for particles i and j. This 
is because [Pij,H] = thus a non-degenerate eigenket 
of H must also be an eigenket of Pij , and the fact (j>o is 
non-negative ensures the eigenvalue can only be -f 1. 

We can now prove that the ground state wave function 
of the spin-1/2 bosonic system must take the form: 

V'^(ri,cri, • • •, Tat, (Tat) = (?!)o(ri, • • ■ , rjv)XAf/2(c^i, • ■ -crjv), 

(3) 

where XAr/2('^ii ' ' ' "'A') is a- spin wave function with quan- 
tum number Stot = N/2. First of all, ipQ is symmetric 
under permutation, since both 0o and Xn/2 are sym- 
metric under permutation. Secondly, tpg is an eigen 
wave function of H with eigen value Eq, thus it must 
be one of the ground states. Lastly but most impor- 
tantly, we need to prove the lowest energy state with 
other Stot — S < N/2 quantum numbers have energies 
higher than Eq . Without losing generality we can choose 
to focus on states with Si^^ = Stot = S. Such a state 
takes the following form: 

"01 (ri: f^i, ■ • • , i"Af, ctat) = ^ (/'5(rp(i), • • ■ , rp(jv)) 

p 

X I Tp(1) • • • Tp(A'i); ip{Ni + l) ■ ■ ■ ip{N)), (4) 

where P stands for permutation, and A^i = N/2 + 5 is 
the number of up spin particles. In order for -0^ to be an 
eigen wave function of H with eigen value E, (pg must be 
an eigen spatial wave function of H with eigen value E, 
which can be chosen to be real. Now we prove E > Eq. 
Since ipg is a common eigen ket of S^^^ and S^^^ with 
S^^f = Stot, we have 

s+Ml) = J2s+\^§)^o. (5) 

i 

Since S^^^ does not annihilate any term of |05 ) expanded 
as in eq. (Q), and matrix elements of S^' are real and pos- 
itive in the S^ representation, in order for eq. (^) to be 
valid 4>s niust contain both positive and negative parts. 
Thus E > Eq. We have thus proved that the ground state 
of the bosonic system described by the Hamiltonian (|l|) 
must have Stot = N/2 with degeneracy 25*40* + 1 = 



if there exists an infinitesimal magnetic field, the degen- 
eracy will be lifted and the spins will be fully-polarized 
along the direction of the field. 

Using similar methods one can also prove the same is 
true for lattice Hamiltonians of the form: 

= - ^ tij{b\„hj„ + h^j^bia) + ^ UiUi + ^ VijUiUj, 

(6) 

where hi„ is the boson annihilation operator for site i and 
spin component cr, tij are real and positive hopping ma- 
trix elements, Vij is the two-body interaction potential, 
and Ui = b\^hi-^ -\- b\ybii is the total boson number on site 
i. The lattice structure is not important and hence not 
specified. 

Wc emphasize our proof does not assume the presence 
of BEC; in fact it is valid in the absence of BEC as well, 
and we will discuss an example of this situation later. 
On the other hand it is crucial that the Hamiltonian H 
is real, so that the eigen wave functions can be chosen to 
be real; the proof is no longer valid when the Hamiltonian 
H is complex, for example when the bosons are charged 
and moving in a magnetic field (with no Zeeman coupling 
so that the SU(2) symmetry is present), or the bosons are 
placed in a rotating trap. 

A physical consequence of the ferromagnetism of the 
ground state is the presence of ferromagnetic spin wave 
excitations. In the following we study the spin wave spec- 
tra of two different situations, one in the continuum and 
the other in a lattice model. 

Consider a continuum system described by the Hamil- 
tonian Eq. with U = 0, so that the system is transla- 
tionally invariant and all eigenkets of H can be chosen to 
have a momentum quantum number k. A natural varia- 
tional wave function for a ferromagnetic spin-wave state 
takes the form 

|k)=5k|0), (7) 
where |0) is the ground state with Stot — Stot — 

N/2, 

and 

S^={l/^)Y,e^'^-^^Sj. (8) 



This is, of course, a straightforward generalization of the 
single-mode-approximation (SMA) for spinless bosons 
Jl2t . Thus a variational approximation (or upper bound) 
for the spin-wave spectrum is 

E{k)=f{k)/s{k)^^, (9) 

where 

/(fc) = (k|(i?-i?o)|k) = ^ (10) 
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is the same oscillator strength as in SMA, while here 

s(fc) = (k|k) = (0|5VklO> = l (11) 

is the spin structure factor. In the SMA for spinless 
bosons the collective mode spectrum is linear for small 
k because the density structure factor is linear in k for 
small k ; here the spin- wave spectrum is quadratic as 
expected for Heisenberg ferromagnets, and we find that 
it is bounded from above by the single particle spectrum. 

For bosons moving in a periodic potential p3[ , the sys- 
tem can be appropriately described by a tight-binding 
lattice Hamiltonian of the form (^). In particular, if the 
periodic potential is strong enough the system loses BEC 
and becomes a Mott insulator Q . In this case the charge 
degree of freedom of the bosons are frozen, but the spin 
degree of freedom are still active and as we show below, 
they can be described by the Heisenberg model with fer- 
romagnetic couplings. For simplicity we consider a boson 
Hubbard model with nearest neighbor hopping and on- 
site repulsion, on a simple cubic lattice: 

<ij>.,a i 
<ij>,a i 

+ '^Un,{n,-l)/2 + const., (12) 

i 

and assume the number of bosons is the same as the num- 
ber of lattice sites, thus there is one boson on every site 
in average. For very large U/t the bosons can not hop 
from one site to another due to the large potential en- 
ergy cost, and the system is in the Mott insulator phase 
with no BEC. A straightforward second-order perturba- 
tion calculation int/U yields an effective Heisenberg spin 
Hamiltonian for the system: 

if = J E (S. • S, + 3/4) (13) 

with ferromagnetic coupling: 

J = -At^/U. (14) 

And the spin-wave spectrum takes the exact form 

E{\i) = {8t^ /U){3 - cos k^- cos ky- cos k^) (15) 

in this limit. This is exactly the opposite for the Hub- 
bard model at half- filling for electrons, where the large 
U/t limit leads to a Heisenberg antiferromagnet. We note 
in passing that in electronic systems, ferromagnetism is 
usually associated with itinerant electrons while localized 
electrons usually give rise to antiferromagnetism; thus 
two-component bosons trapped to lattice potentials pro- 
vide new opportunities to study lattice ferromagnets. 



We now turn our discussion to the relation between our 
work and existing theoretical work on this and related 
fields. 

Ho and Yin |^ considered the general spin structure of 
Bose-Einstein condensates of bosons with arbitrary spin. 
In their work they assumed weakly interacting bosons 
thus to a very good approximation, all bosons occupy 
the lowest-energy orbital state. They pointed out for 
spin- 1/2 bosons due to Bose statistics, the bosons must 
also occupy the same spin state and the system is a 
ferromagnet, and very appropriately termed the ferro- 
magnetic state "statistical ferromagnet". Our work has 
considerable overlap with theirs in spirit, and certainly 
comes to the same conclusion as theirs. On the other 
hand our proof is rigorous and more general; it applies 
to interacting bosons with arbitrary interaction strength 
and in particular, strongly interacting bosonic systems 
without BEC. We have thus generalized the statistical 
ferromagnet to a much wider class of bosonic systems. 

Rojo studied two-component bosonic systems with 
component- dependent interactions, in which the inter- 
component interaction is always more repulsive or less at- 
tractive than the intracomponent interaction. He showed 
that in this case the ground state is "fully polarized" in 
the sense that only bosons of one component appear. 
In the spin analogy the interaction he considered has 
Ising anisotropy, and the system is an Ising ferromagnet. 
Thus his work is complementary to ours, which studies 
a Heisenberg ferromagnet. By the same token, in his 
model if the intercomponent interaction were less repul- 
sive or more attractive than the intracomponent interac- 
tion, the interaction would have XY anisotropy and the 
system would be an XY ferromagnet. He further found 
that as the two interactions become closer and closer to 
each other (or one is approaching the Heisenberg point), 
the velocity of one of the two linear collective modes van- 
ishes. This, as we see here, is because at the Heisenberg 
point the spin wave spectrum has quadratic dispersion at 
long wave length. 

Li and coworkers have studied one-dimensional 
SU(2) bosons with (5-function interactions. They used 
Bethe ansatz to obtain the ground and low-lying state 
spectra and quantum numbers, and showed the ground 
state is fully spin polarized and there is a quadratic spin 
wave mode. Their results certainly agree with the general 
conclusions of the present work. 

Pu et al. 1^ studied possibility of ferromagnetic states 
in bosonic systems originated from dipole-dipole inter- 
actions. In their model the ferromagnetism comes from 
interactions that are magnetic in the first place, which is 
very different from our case, where the ferromagnetism 
is essentially kinetic energy driven. 

Thus far, our discussions have been focusing on bosonic 
systems. Fermionic systems are of course very different, 
due to the difference in statistics. On the other hand 
in two-dimensions (2D) the statistics of the particles can 
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be transformed through flux attachment transformations 
psf . Thus our results also have implications on elec- 
trons in a strong magnetic field in the quantum Hall 
regime. For example, it gives a new perspective to the 
spin/pseudospin ferromagnetism for single/bi-layer quan- 
tum Hall systems at Landau level filling factor aX v = 1 
[^,0, as the electrons can be mapped to bosons carry- 
ing one flux quantum, and the flux carried by the bosons 
cancel that of the external magnetic field in average, thus 
at the mean-field level these systems are mapped onto bo- 
son systems at zero magnetic field, precisely the system 
we study here! Using the same transformation one can 
also map spin- 1/2 bosons in a magnetic field with u ^ 1 
to fermions with zero field, where one expects a singlet 
ground state. This has indeed been found to be the case 
in finite-size studies, for Coulomb interaction [|8|. 
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